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ABSTRACT 

Aims. HESS J1858+020 is a weak 7— ray source that does not have any clear cataloged 
counterpart at any wavelengths. Recently, the source G35.6-0.4 was re-identified as a 
SNR. The HESS source lies towards the southern border of this remnant. The purpose 
of this work is to investigate the interstellar medium around the mentioned sources in 
order to look for possible counterparts of the very-high energy emission. 
Methods. Using the 13 CO J=l-0 line from the Galactic Ring Survey and mid-IR data 
from GLIMPSE we analyze the environs of HESS J1858+020 and SNR G35.6-0.4. 
Results. The 13 CO data show the presence of a molecular cloud towards the southern 
border of SNR G35.6-0.4 and at the same distance as the remnant. This cloud is 
composed by two molecular clumps, one, over the SNR shell and the other located 
at the center of HESS J1858+020. We estimate a molecular mass and a density of 
~ 5 x 10 3 Mq and ~ 500 cm -3 , respectively for each clump. Considering the gamma- 
ray flux observed towards HESS J1858+020, we estimate that a molecular cloud with a 
density of at least 150 cm -3 could explain the very-high energy emission hadronically. 
Thus, we suggest that the 7— ray emission detected in HESS J1858+020 is due to 
hadronic mechanism. Additionally, analyzing mid-IR emission, we find that the region is 
active in star formation, which could be considered as an alternative or complementary 
possibility to explain the very-high energy emission. 
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1. Introduction 

Among the rich population of TcV 7— ray sources, a good fraction are associated with 
Galactic phenomena which include supernova remnants (SNRs), pulsar wind nebulae 
(PWNc) and binary systems (Hinton 2008). Extragalactic TeV sources are associated with 
active AGNs (mostly blazars). There are, however, a number of 7— ray sources for which 
there are not yet a clear identification and an origin established. In this paper we focus on 
the very-high energy source HESS J1858+020. 
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HESS J1858+020 is a weak 7-ray source detected with H.E.S.S. (the High Energy 
Stereoscopic System), an array of air Cherenkov telescopes. Though nearly point-like 
source, its morphology shows a slight extension of ~ 5' along its major axis. The source 
has been detected at a significance level of la with a differential spectral index of 2.2 
±0.1. At the present, HESS J1858+020 docs not have any clear cataloged counterpart at 
any wavelengths. The pulsar PSR J1858+0143, which is energetic enough to power HESS 
J1858+020, is located far away from the center of the gamma source, thus an association 
between them is unlikely (Aharonian et al. 2008a). 

Recently Green (2009), using radio continuum data from the Very Large Array (VLA) 
Galactic Plane Survey (VGPS; Stil et al. 2006), re-identified the radio continuum source 
G35.6-0.4 as a supernova remnant (SNR). This remnant is seen in projection over the 
northern border of HESS J1858+020. Green (2009) estimated an age of 30000 years old 
for the SNR and a distance of ~ 10.5 kpc based on the proximity of the remnant with 
the HII region G35. 5-0.0. As pointed out in several works (see Aharonian & Atoyan 1996; 
Yamazaki et al. 2006; Gabici et al. 2007, 2009), a SNR, in particular an old remnant (i.e 
age larger than a few 10000 years), interacting with a molecular cloud could explain the 
origin of the 7— ray emission via pion decay from proton-proton collisions. In this context, 
we expect a correlation between 7— ray emission and matter concentration. 

In this letter we study the interstellar medium around the sources G35.6-0.4 and HESS 
J1858+020 to investigate the possible origin for the very-high energy emission. 

2. Data 

To analyze the ISM towards HESS J1858+020 we used the 13 CO J=l-0 emission obtained 
from the Galactic Ring Survey. This survey maps the Galactic Ring with an angular and 
spectral resolution of 46" and 0.2 km s _1 , respectively (see Jackson et al. 2006). We also 
used the mosaiced images from GLIMPSE and MIPS GAL and the GLIMPSE Point-Source 
Catalog (GPSC) in the Spitzer-IRAC (3.6, 4.5, 5.8 and 8 ^m). IRAC has an angular 
resolution between 1'.'5 and L"9 (see Fazio et al. 2004 and Werner et al. 2004). MIPSGAL 
is a survey of the same region as GLIMPSE, using MIPS instrument (24 and 70 /zm) on 
Spitzer. The MIPSGAL resolution at 24 fj,m is 6". 

3. Results and Discussion 

Analyzing the whole 13 CO J=l-0 data cube in a region containing SNR G25.6-0.4 and 
HESS J1858+020, we find a molecular cloud lying in the southern border of the remnant 
in the velocity range between 51 and 59 km s . This velocity range coincides with the 
systemic velocity of G25.6-0.4 derived from its estimated distance (Green 2009). Figure 1 
shows a three color image where red represents the Spitzer-IRAC 8 /im band, green is the 
MIPSGAL emission at 24 /im and blue is the radio continuum emission of G35.6-0.4 at 
20 cm from the VGPS. Black contours delineate this last emission and white contours are 
the 13 CO ,1=1-0 emission integrated between 51 and 59 km s" 1 . As mentioned above, the 
extension of the HESS source is slightly elliptical with ~5' along its major axis. As the 
extent and angle are not well constrained, in Figure 1 we represent the position and the 
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extension of HESS J1858+020 with a circle of 5' in diameter. From this figure, it can be 
seen that the molecular cloud presents two well defined clumps, one, centered at I = 35?60, 
b = -0?53, over the SNR shell, and the other, centered at I = 35?58, b = -0?58, in 
coincidence with the center of HESS J1858+020. The 7— ray emission could be related to 
both molecular clumps. 

By inspecting the 13 CO spectra we found that the clump located over the SNR shell 
shows some possible kinematical evidence of shocked gas. As Fig. 2 displays, the spectra 
are not symmetric and they present a slight spectral shoulder or a less intense component 
at "blucshiftcd" velocities. It could be evidence of turbulent motion in the gas, may be 
produced by the SNR shock (see e.g. Falgarone et al. 1994). Taking into account the 
positional and kinematical agreement between the molecular gas and the SNR G35.6-0.4, 
we suggest that the remnant is interacting with the adjacent molecular cloud. 




HESS J1858+020 
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Fig. 1. Three color image of the analyzed region (red = 8 /im, green = 24 /im, and blue = 20 
cm radio continuum emission). The black contours delineate the radio continuum emission with 
levels of 20 and 28 K. The white contours are the 13 CO J=l-0 emission integrated between 51 
and 59 km s , its levels are 3.2 and 5 K km s _1 . The yellow circle represents the source HESS 
J1858+020. 



Using the 13 CO J=l-0 line and assuming local thermodynamic equilibrium (LTE) we 
estimate the H2 column density towards the molecular clumps described above. We use: 



T cx fr 13 dv 



N( i3 CO) = 2.42 x nr 



1 - exp(-5.29/T cx ) 

to obtain the 13 CO column density. T13 is the optical depth of the line, we assume that T ex is 
10 K and the 13 CO emission is optically thin. We use the relation N(H 2 )/N( 13 CO)~ 5 x 10 5 
(e.g. Simon et al. 2001). Finally, we estimate a molecular mass and a density of ~ 5 x 10 3 
M.q and ~ 500 cm -3 , respectively, for each molecular clump. The mass value was obtained 
from: 

M = / im H ^[fl 2 flN(H 2 )], 
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Fig. 2. 13 CO J=l-0 spectra obtained towards the maximum of the molecular clump that lies over 
the SNR shell. The spectra, which show brightness temperature vs. velocity, correspond to nine 
positions observed in a region of about 50" x 50" centered at the molecular clump peak. 

where il is the solid angle subtended by the 13 CO J=l-0 beam size, run is the hydrogen 
mass, fi, the mean molecular weight, is assumed to be 2.8 by taking into account a relative 
helium abundance of 25 %, and D is the distance assumed to be 10.5 kpc. Summation was 
performed over all the observed positions within the 3.5 K km s _1 contour level (not shown 
in Figures 1 and 3). 

We can estimate the required density matter in the 7— ray production region for hadrons 
from the observed 7— ray flux of HESS J1858+020. From 

<*F 7 /dE(> E mm ) = JV (£/lTcv)- r , 

we obtain ~ 1.2 x 1(T 12 cm^s" 1 . Here E mm = 0.5 TcV, V = 2.17, N = 0.6 x 10~ 12 
cm _2 s _1 (Aharonian et al. 2008a). Using the equation 16 in Torres et al. (2003), we obtain a 
density of about 150 cm~ 3 assuming an acceleration efficiency of hadrons of the order of 3% 
and a supernova power of 10 51 ergs. Thus the molecular gas in positionally coincidence with 
HESS J 1858+020 is densely enough to generate the very- high energy emission hadronically. 

Regarding the IR emission in the studied region, from Fig. 1 it can be seen the presence 
of at least two bright sources embedded in the southern molecular clump. Also the mid-IR 
emission shows the presence of hot dust and policyclic aromatic hydrocarbons (PAHs), 
green and red in Fig. 1, respectively, suggesting an active star forming region. In order 
to confirm that, we carried out an IR photometric study of the sources that lie over the 
molecular clump that is centered at the HESS source position. 

In Figure 3 we present a three color image with three IRAC-Spitzer bands: red is 8 
jum, green is 4.5 /im and blue is 3.6 /^m. As in Figure 1 the white contours are the 13 CO 
J=l-0 emission integrated between 51 and 59 km s _1 . To look for tracers of star formation 
activity, we use the GLIMPSE Point-Source Catalog to perform photometry. Considering 
only sources that have been detected in the four Spitzer-IRAC bands, we found 32 sources 
within a circular area of 80" in radius centered at the molecular clump. According to 
the criteria presented by Allen et al. (2004), we found 6 sources as young stellar objects 
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Fig. 3. Three color image of IRAC-Spitzer bands (red = 8 /im, green = 4.5 pm, and blue = 3.6 
/im). The white contours are the 13 CO J=l-0 emission integrated between 51 and 59 km s , 
its levels are 3.2 and 5 K km s" 1 . The crosses are YSO candidates according to the photometric 
study described in the text. Green crosses are Class I sources and cyan crosses represents Class II 
sources. 



(YSO) candidates. Two of them are Class I objects (protostar with circumstellar envelope; 
green crosses in Fig. 3) and the others are Class II objects (young stars with only disk 
emission; cyan crosses in Fig. 3). From Fig. 3 it can be noticed a bright source that lies 
almost at the center of the molecular clump, but in this case since it has some bands with 
null fluxes in the GLIMPSE catalog, we can not perform photometry. However, from Fig. 
3 can be appreciated that this source appears slightly extended in the 4.5 /im emission 
(green), which according to Cyganowski et al. (2008), may suggest that it could be a 
MYSO driving outflows. For the sources that have also detection in the 2MASS JHK- 
bands, we derive their spectral energy distribution (SED) by fitting the fluxes using the 
tool developed by Robitaille et al. (2007) and available online 1 . This was possible only 
for three Class II sources, resulting that they are indeed young objects (~ 10 5 years). 
We conclude that this region is active in star formation, which can be considered as an 
alternative or complementary possibility in order to explain the very-high energy emission. 
HESS J 1858+020 could be a similar case as SNR W28. Several molecular clumps and 
star-forming regions were found towards W28 and proposed to be candidates to explain 
the very- high energy emission detected in the region (Aharonian et al. 2008b). Further 
sub-mm observations, high CO transitions and tracers of star formation, will provide more 
accurate mass and density estimates and allow to search for perturbed gas by both SNR 
shock and star formation processes. 



http:/ /caravan. astro.wisc.edu/protostars/ 
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4. Summary 

In this work we investigate the ISM towards the southern border of the recently re-identified 
source G35.6-0.4 as a SNR, which coincides with the northern portion of the very-high 
energy source HESS J1858+020. Analyzing the 13 CO J=l-0 emission we find a molecular 
cloud in positional and kincmatical agreement with the SNR, suggesting that the molecular 
gas is being affected by the remnant shock. The discovered molecular cloud is composed 
by two clumps, each with mass and density of ~ 5 x 10 3 M© and ~ 500 cm -3 , respectively. 
Considering the gamma-ray flux observed towards HESS J1858+020, we estimate that a 
molecular cloud with a density of at least 150 cm~ 3 could explain the very-high energy 
emission hadronically. Thus, we suggest that the interaction between the SNR G35.6-0.4 
and the molecular cloud is responsible for the 7— ray emission. Additionally, analyzing mid- 
IR emission, we find that the region is active in star formation, which could be considered 
as an alternative or complementary possibility to explain the very-high energy emission. 
This SNR - molecular cloud - HESS source complex could be a similar case as the high- 
energy sources detected towards SNR W28. Further sub-mm observations will provide more 
accurate mass and density estimates and allow to search for perturbed gas by both SNR 
shock and star formation processes. 
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